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ABSTRACT 

The determination of the critical exponents by means of the Exact Renormalizion 
Group approach is still a topic of debate. The general flow equation is by con- 
struction scheme independent, but the use of the truncated derivative expansion 
generates a model dependence in the determination of the universal quantities. 
We derive new nonperturbative flow equations for the one-component, Z2 sym- 
metric scalar field to the next-to-leading order of the derivative expansion by 
means of a class of proper time regulators. The critical exponents 77, v and u 
for the Wilson-Fisher fixed point are computed by numerical integration of the 
flow equations, without resorting to polynomial truncations. We show that by 
reducing the width of the cut-off employed, the critical exponents become rapidly 
insensitive to the cut-off width and their values are in good agreement with the 
results of entirely different approaches. 
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Universality is the property which guarantees that, at a critical point, the di- 
vergences of the correlation functions of an entire class of physical systems have 
the same power law behavior with the same exponents (critical exponents), re- 
gardless of their specific features and especially of their microscopic (ultraviolet 
(UV)) details. This notion is very general and it often appears in the study of 
critical phenomena and in quantum field theory. The prediction of those expo- 
nents has been recognized to be one of the most important achievements of the 
renormalization group approach pioneered by Wilson and others [0], 0]. A stan- 
dard approach to this problem is the e-expansion [[| method, whose main feature 
is an expansion in e = 4 — D (where D indicates the number of dimensions) of the 
renormalized free energy. Unfortunately this method is not entirely satisfactory, 
since actual values for three-dimensional systems are obtained only after setting 
e = 1 at the end of the calculation. On the other hand, ordinary perturbation 
theory in such nonperturbative framework requires quite involved loop calcula- 
tions ||. Moreover since both e-expansion and perturbative expansion at fixed 
dimension yield divergent series, one has to rely on the Borel summability of these 
series in order to get meaningful results. 

In recent years the Exact Renormalization Group (ERG) approach has proven 
to be very useful in many situations where the usual perturbative methods are 
not totally satisfactory [^j, [5]. The main step in this approach is to solve the flow 
equation which is a complicate non-linear functional equation. For this reason one 
must then resort to sofisticate approximation methods in order to extract physical 
information. The derivative expansion allows a systematic expansion of the renor- 
malized action in terms of local field operators. It has the advantage of recasting 
the functional flow equation into a set of coupled partial non-linear differential 
equations for the local potential £4(0) and for the wavefunction renormalization 
function Zk(4>) in the next-to- leading order approximation which improves the so 
called Local Potential Approximation (LPA) where only the equation for Uk(4>) 
is considered and Zk(<p) = 1 is kept fixed. This topic has been the subject of 
an intense investigation and the reader is referred to the excellent review in || 
for a recent survey of results for the scalar theory. Although, as far as the cal- 
culation of universal quantities is concerned, the exact evolution equation is by 
construction scheme independent, this no longer holds in specific truncations of 
the full evolution equation. In fact different realizations of the smooth cut-off 
ERG equations determine critical exponents which are different if, for instance, 
an exponential cut-off or a power law cut-off is used || |9|. 

Recently, the ERG scheme dependence has been the subject of many papers. 
In |H| the scheme dependence has been analyzed in order to optimize the con- 
vergence of approximate solutions of the ERG equations. Other authors have 
considered the effects of the scheme dependence in the convergence of polynomial 
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truncations of the flow equation in the LPA [O, IT2[. Actually it would be de 



sirable to have an approximation framework which allows to control the scheme 
dependence to each order in the derivative expansion, without resorting to further 
polynomial truncations and which, at the same time, could provide an estimate 
of the uncertainty on the determination of the universal quantities. 

The first determination of the critical exponents u, u from a numerical reso- 
lution of the flow equation has been performed in |13j by using the LPA of the 
Wegner-Houghton (WH) equation for the wilsonian action, [14], which is a sharp 
cut-off version of the ERG equation. However the WH equation, truncated to 
the next-to-leading order in the derivative expansion although predicts the cor- 
rect two loop value of the field anomalous dimension 7] in D = 4 |]15| , fails to give 
a reasonable value of i] at the Wilson-Fisher fixed point in D 



3 16 



To overcome the difficulties with the sharp cut-off formulation, several authors 



17] , |T8| , |H| have proposed a smooth cut-off realization of the ERG equations by 
means of the Schwinger proper time regulator. This approach has been used 
in many contexts like Yang-Mills theories |2(| and Chiral Symmetry Breaking 
In |12[ the effect of this regulator on specific polynomial truncations of 

H 
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the LPA has been discussed. In 



the same cut-off function is used to include 
the effect of a field independent wavefunction renormalization function Zk in the 
O(N) model. However the inclusion of the field dependence in Zk turns out to be 
essential for a consistent determination of the scaling properties of the system. 
A discussion on the determination of the wavefunction renormalization function 
within this framework has been presented even in [TJ3|], but so far the full flow 
equation for Zk(4>) within this approach has not been derived. 

The aim of this paper is to derive the full coupled flow equations for Z^cj)) 
and Uk(4>) by means of the general class of proper time regulators introduced in 



17] , [I8| , |19|| and to numerically determine the anomalous dimension rj and the 



critical exponents v and to, by solving these equations around the Wilson-Fisher 
fixed point in D = 3. We also show that, when diminishing the size of the cut-off 
width, the critical exponents become stable, converging to fixed values. 

Let us now review the basic assumptions in the derivation of our evolution 
equation. The one-loop contribution to the effective action is given by the trace 
of \og(5 2 S/ (5(j)5(j))), where S is the classical action. By introducing the Schwinger 
proper time representation for the logarithm with the proper time variable s, we 
can express the one-loop contribution as 



q1 — loop 
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where we have introduced the heat-kernel smooth regulator f^. The usual one- 
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loop effective action is recovered by setting f k = 1 and properly regulating the 
UV divergences in Eq. flU). 

In the spirit of the derivation of the WH equation, the flow of the wilsonian 
action is obtained by taking as a sharp momentum cut-off, i.e. a Heaviside 
function which is equal to one for momenta greater than k and vanishes for 
momenta smaller than k, and then by performing the integration of the modes 
of the action through infinitesimal steps 5k in the momentum space. Therefore 
the action Sk-sk comes from the integration of the modes between k — 5k and k 
in the action defined at the scale k. A smooth cut-off on the momenta in this 
framework is then obtained by choosing a non-singular instead of the Heaviside 
function and the differential flow equation reads 

dS k 1 r°° ds ,df k i 5 2 S k 



1 [°° ds df k , 5 2 S k ^ 
"2 Tr Jo T^ eXP i~W ( ) 



k dk 

It must be remarked that the WH equation contains terms that are neglected 



in Eq. (E3) but, as discussed in [23], these terms do not contribute to the next-to- 



leading truncation in the derivative expansion. 



About f k one can in general require the following properties (see [18|, (jlj): 1) 
f k must depend on k through the dimensionless quantity sk 2 . We shall consider, 
as in p2]| the variable sZ k k 2 , where Z k is the wave function renormalization 
function, in order to reproduce the usual form of the propagator in the following 
equations. 2) f k=0 = 1 so that S k=0 must be identified with the full effective 
action. This implies that in the limit k — > the infrared cut-off is removed. 3) 
f k (sZ k k 2 ) must tend to zero sufficiently rapidly for large sZ k k 2 to suppress the 
small momentum modes. In particular we choose, for any integer m > 1 

ft\sZ k k 2 ) = e~^ k2 Y, { -^T 1 (3) 

i=0 % - 

We note that, as it stands, Eq. (0) for some values of m and of the dimension 
D is not UV convergent, but in D = 3 and for m > 1, which is the case we are 
interested in, no UV divergence appear in the differential flow equation. 

An analysis of the polynomial truncations of the LPA to Eq. (0) for m — 1, 2, 3 
has been performed in [|l^], where it is shown that all the parameters extracted 
from the flow equations, regarded as functions of the polynomial truncation, 
converge more rapidly when m grows. This is a first indication that the RG 
equations provide more stable results for increasing values of m. 

As anticipated above, we shall focus here on the first order truncation in 
the derivative expansion of the action S k which corresponds to the following 
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parametrization 

S k = J d D x(^-d^dy + V k {<t>)) (4) 

and therefore Eq. (g) is reduced to two coupled partial differential equations for 
Vk(4>) and Zk((j)). No further polynomial expansion or truncation in powers of 
the field is made on Vk(4>) and Z k ((f)). 

In order to project the flow equation on the local potential term and on the 
kinetic term we first set (see also |24], [23], [15)) <fi = <po + (p(x) where 0o is a constant 



field configuration and f(x) is a small fluctuation field. Next, we expand V k (4>) 
and Z k (<p) in both sides of Eq. (fj) and retain only terms at most quadratic in 
dip. Then we define 5 2 Sk/ (5<j)5<j)) = A = Aq + 5 A where Aq is only a function of 
the constant field <po and 5 A contains all the linear and quadratic contributions 
in dip. Thus the trace becomes 

_I Tr r ds k dh e _ sA = 

2 Jo s ok 

I d D x r ds Z k k 2 (sZkk p m e- sZ ^(x\e-^ Ao+5A ^\x) = 
J Jo ml 



ds Z h k A k . ] e- sZkk2 (x\e- sAo (l - s5A 
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+^{[5A,A ]+5A 2 ]} + ...)\x) (5) 

where in order to disentangle the trace in Eq. (|5|) we used the Baker-Campbell- 
Hausdorff formula, and the dots stand for the higher order terms in the s expan- 
sion of the exponential. Actually, in order to collect all the contributions to the 
wavefunction renormalization function we had to retain terms up to the 0(s 4 ) 
order in this expansion. The evaluation of the trace, as explained in fl24| , ^3| , |i~5|j , 
is then performed by inserting the identity / dp \p)(p\ in the trace in Eq. (|^) 
and by making use of the commutation relation [p M ,/(x)] = id fl f(x). Finally, 
by collecting on both sides of Eq. ([5]) the coordinate independent terms and the 
coefficients of the kinetic term dipdip, we get respectively the flow equation for 
and for Z^. Remarkably, all the integrals in the proper time s can be performed 
analytically and therefore we end up not with integrodifferential flow equations, 
but simply with differential equations. 

It is convenient to express the flow equations in terms of new dimension- 
less quantities t, V(t, $), Z(t,&), defined by the following relations: t = 
ln(A/fc), (j) = k^l 2 ^a^ } V k (k,<j)) = ak 3 V(t,®) and Z k (k,<j>) = fc-"Z(t, $) 
where rj is the anomalous dimension of 0, A is a UV reference momentum scale, 
and the constant a is defined in terms of the Gamma functions: a = (T(m — 
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l/2))/(87r v /7rm!)). The flow equations then read (the / indicates derivative w.r.t. 
the field $) 

^ = 3V--i±^'-pU m - 1/2 (6) 
dt 2 \Z + V"J V ; 

l±l*x + f-g-) "" 1/2 f (m - 1/2) (z« - « 

' 2 ^Z + 1/'7 V Z + V" V 24Z 

7K - 1/4) 7VW , (m 2 -l/4)(m + 3/2) 2 \ 
6(Z + V") 2 Z V + 6(Z + V"f Z[V ] ) (?) 

In order to determine the critical exponents we need to find the fixed points 
of Eqs. (|||7|) and then to study the linearized version of the equations around the 
fixed point solutions. The latter are the t-independent solutions of Eqs. (|6],|7|). 
They are indicated as V*(<&) and Z*($). Once V*(<&) and Z*(<3>) are determined, 
one can linearize Eqs. (§,0), by introducing small t-dependent perturbations 
around V*(<&) and Z*($): 

v(t, $) = v*($) + $v(t, $) = + e A '^($) 

Z(t, $) = Z*($) + $) = Z*($) + e xt z($>) (8) 

and by retaining in Eqs. on ly linear terms in v(<&), z(§). 

It is worthwhile to notice that Eqs. (§,0) are invariant under the transfor- 
mations $ — > e$, Z — > e~ 2 Z. This is the consequence of a reparametrization 
invariance of Eq. (|2|) where and Zk transform as $ and Z respectively and the 
proper time is redefined as s — > e 2 s, and all the other quantities are not modified. 

Before going on we consider in more detail the behavior of the regulator /& 
when m grows. As a simple check one can go back to Eq. (|5|) and neglect the Z^ 
effects. In this case 5 A = and no expansion of the exponential is needed and 
only the term exp(— sAq) is left. By integrating the variable s we get 
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Clearly the p integration leads back to Eq. (Q) with Z = 1, but from the 
integrand in the r.h.s. of Eq. (||) it is possible to learn which momenta p are 
relevant in the integration. In Fig. 1 the rescaled integrand F(p) = mp 2 (k 2 j (k 2 + 
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p 2 + Vk")) m+1 for some values of m and for fixed k is plotted. The maximum 
of F(p) is at p 2 = (k 2 + V k ")/m and F(p) = (k 2 + V k ")- m (mk 2 /(m + l)) m+1 . 
So, when m grows, the cut-off f k has the effect of selecting smaller and smaller 
momentum shells, centered at p where the integrand behaves as a power of the 
propagator. Then, the large m limit is a kind of sharp cut-off limit for a modified 
propagator, which however cannot be directly related to the sharp cut-off WH 
equation. Indeed the logarithmic behavior of the latter is clearly different from 
the one of Eq. (|6|) where the term (Z/(Z + V")) rn ~ 1 / 2 appears instead. At the 
same time no straightforward limit m — > oo can be taken in Eqs. (§,0) in order 
to determine the analytic asymptotic behavior of the critical exponents. This can 
be seen in the LPA from Eq. ([]) with Z = 1. In this case, when m — > oo, the 
sign of V" determines whether the last term in (|6|) diverges or vanishes (or it is 
one for V" = 0). Therefore no uniform convergence of the solution is found in 
the limit m — > oo. Unfortunately, a deeper understanding of the relation between 
this heat-kernel regulator and the other formulations of the ERG equations is still 
missing and we have to limit ourselves to compare the numerical results of the 
different methods. 

We are interested in determining the anomalous dimension rj at the only non- 
gaussian fixed point solution, the Wilson-Fisher fixed point, and the exponents 
v and u which are related to the eigenvalues A of the linearized equations, being 
defined as the inverse of the only positive eigenvalue and as the opposite of the 
less negative eigenvalue respectively. 

The first step concerns the determination of non-trivial solutions V*($) and 
Z*($) and we closely follow the procedure outlined in |9j (see also [[RJ). We 



use the shooting method in order to numerically solve the fixed point equations 
taking, as boundary conditions, the asymptotic expression of the potential and 
wave function renormalization function at large $ together with the symmetry 
property of these functions at the origin: Z*'(0) = V*'(0) = 0. The normalization 
at the origin Z*(0) must be fixed to some definite value and, as a test of the 
reparametrization invariance, we checked that the numerical results found for the 
critical exponents do not depend on the particular value of this normalization. 
The procedure to determine the asymptotic expressions from the equations is 
explained in detail in || and therefore we do not repeat it here. In practice, 



as it happens for the sharp cut-off equations, ||16|| , we find that the numerical 
resolution of the problem is simpler if the equation for the potential is replaced 
by the corresponding equation for its field derivative V*'($>). For each value of m 
we find only one non-trivial fixed point as it was expected and the values of the 
anomalous dimension for each m are displayed in Table 1. 

Having determined the fixed point, we again use the shooting procedure to 
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m 


V 


m 


V 


1 


0.0653 


11 


0.0365 


2 


0.0507 


12 


0.0362 


3 


0.0452 


13 


0.0360 


4 


0.0423 


14 


0.0358 


5 


0.0405 


15 


0.0356 


6 


0.0393 


16 


0.0354 


7 


0.0385 


17 


0.0353 


8 


0.0378 


20 


0.0350 


9 


0.0373 


30 


0.0343 


10 


0.0369 


40 


0.0340 



Table 1: The anomalous dimension rj at the Wilson-Fisher fixed point determined 
for various values of m. 

evaluate v and u from the linearized equations as explained in ||. The results are 
reported in Table 2 as the 0(d 2 ) approximation estimates. In order to test the 
improvement due to the inclusion of the wave function renormalization function, 
we also show in Table 2 the values of is, uj obtained within the LPA ( O(<9 ) ), i.e. 
from the linearized potential equation with rj = and Z = 1 fixed. 

We note that the numerical procedure employed to solve the equations be- 
comes more problematical and converges less rapidly for larger values of m. More- 
over, since the fixed point equations are more stable than the linearized ones we 
can push m up to m = 40 in the former case but only up to m = 9 in the latter 
case. 



m 


is 


UJ 




O(d ) 


0(d 2 ) 


o(d°) 


0(d 2 ) 


1 


0.6604 


0.6348 


0.629 


0.847 


2 


0.6439 


0.6331 


0.674 


0.762 


3 


0.6381 


0.6311 


0.698 


0.738 


4 


0.6351 


0.6299 


0.711 


0.727 


5 


0.6333 


0.6290 


0.720 


0.721 


6 


0.6322 


0.6286 


0.727 


0.717 


7 


0.6313 


0.6279 


0.731 


0.713 


8 


0.6306 


0.6277 


0.735 


0.711 


9 


0.6301 


0.6276 


0.738 


0.710 



Table 2: Values of the exponents is,u for various m. The 0(d°) (LPA) values 
are obtained from the potential equation only with Z = 1, rj = fixed. 
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The first thing that should be noted is that the change of rj, v, 00 diminishes 
when m grows and the value of these parameters becomes almost constant : v and 
w (0(<9 2 )) already at m = 9 are almost stable up to the third significant digit and, 
according to the trend shown up to m = 40, a rapid convergence of 77 is observed. 
The dependence on the parameter m is practically vanishing. Since the analytical 
approach to the m — > 00 limit of Eqs. (|6],0) is troublesome, we have to rely on 
the numerical results shown in Tables 1,2 in order to understand the behavior of 
77, u, uj for large values of m. As an additional check we have considered a fit to 
the numbers in Table 1 with various analytical trial forms of the function rj(m), 
containing a few free parameters. In these checks the x _s< l uare analysis clearly 
indicates a finite asymptotic value of 77(777. — > 00) around 77 = 0.0329. 

As far as the derivative expansion improvement is concerned, v and u behave 
differently in the two approximations here considered: the former is always de- 
creasing with m; the latter has a different trend in the two approximations and 
the LPA gives a better estimate than the 0(d 2 ) approximation. However it must 
be remarked that both u, uj become more stable with m when both equations for 
v and z are considered, rather than in the LPA. It is conceivable to expect that, 
when going to higher orders in the derivative expansion, the critical exponents 
converge to a value that sits between the 0(d°) and the 0(d 2 ) estimates. 

In Table 3 we report some estimates of the critical exponents obtained by 
various formulations of the ERG equations and by completely different methods 
(see captions for references). Table 3 collects only few relevant results and it is by 
far an incomplete list. Much more exhaustive lists can be found in [RL H, 0, B3] . 
A comparison with the results in Tables 1 and 2 indicates that the particular 
cut-off here employed provides a little improvement in the agreement with the 
non-ERG approaches, on the anomalous dimension determination. Concerning v 
and u) our determinations are within the average, but slightly smaller (although 
it should be noticed the small experimental value of v) . 

Therefore, the convergence of the numerical results in Tables 1 and 2 to such 
good values in the limit of large m, provides an indication that the heat-kernel 
cut-off, in that limit, becomes a sensible regulator of both ultraviolet and infrared 
modes. Moreover it is particularly effective in reducing the weight of all the 
irrelevant operators neglected in our derivative expansion truncation, which, on 



the other hand, as it is argued in |nj , become important for the same truncation 
of the WH equation. 

In conclusion, we have presented a particular version of the RG equations for 
a single field scalar theory with a heat-kernel cut-off, introduced in the Schwinger 
proper time formalism. This cut-off is parametrized by an index that regulates 
the size of the momentum shell integrated in the blocking procedure. We con- 
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V 


V 


us 


(a) 


0.0467 


0.6307 




\°) 


U.UOoy 


U.Olol 


u.oy t 




U.U4U^ ( ) 


u.ozo^y ) 


U.ooy 1 ) 


(r\\ 


n n/io 

U.U4/ 


n rod 


U. 1 04 


(e) 


0.0335(25) 


0.6304(13) 


0.799(11) 


(f) 


0.0360(50) 


0.6290(25) 


0.814(18) 






n conn/ 1 \ 
U.ooUU^loJ 


U.ozo^oU j 


\ n ) 




U.UO 1U yO J 




(i) 


0.026(3) 


0.624(2) 


0.82(3) 


(1) 


0.0374(14) 


0.6294(9) 


0.87(10) 


(P) 




0.625(6) 


0.80(5) 



Table 3: Values of the exponents 7], u, uj, calculated with various methods. 

(a) RG effective average action approach with exponential smooth cut-off, to 
0(d 2 ) order 

(b) ERGE with powerlike smooth cut-off to 0(d 2 ) order f^|/. 



(c) ERGE expansion into a truncated set of ordinary differential equation [26], 
quoted from j^j 



(d) Polchinski version of the ERGE [j27j . 

(e) Seven-loop perturbation series in D = 3 j^j. 

(f) Five-loop e-expansion ^J. 

(g) high temperature series quoted from j^j. 

(h) high temperature series %2Q] , quoted from j^J. 

(i) Monte Carlo simulation fiSQJ . 

(I) Monte Carlo simulation fftD^, quoted from j^J. 

(p) experimental data from the liquid-vapour transition, quoted from jffl/. 



sidered a derivative expansion of the flow equation truncated to the 0(d 2 ) order, 
obtaining a set of two coupled equations which are reparametrization invariant. 
The corresponding determination of i], u, uj at the Wilson-Fisher fixed point in 
three dimensions is substantially encouraging due to the good agreement with 
the other estimates of these parameters. 



We would like to acknowledge Martin Reuter for very enlightening discussions 
and encouragements. 
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Figure 1: The function F(p) (see text), plotted vs. p, for fixed k = 1, V" = 0.01 
and m = 1,2, 10,20. 
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